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ABSTRACT 



We present a simple heuristic model for the time-averaged soft X-ray temperature distribution in the accretion spot on the white 
dwarf in polars. The model is based on the analysis of the Chandra LETG spectrum of the prototype polar AM Her and involves 
an exponential distribution of the emitting area vs. blackbody temperature a(T) = ag exp(-T /T ). With one free parameter besides 
the normalization, it is mathematically as simple as the single blackbody, but is physically more plausible and fits the soft X-ray and 
far-ultraviolet spectral fluxes much better. The model yields more reliable values of the wavelength-integrated flux of the soft X-ray 
component and the implied accretion rate than reported previously. 
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1. Introduction 

In the 'polar' subtype of cataclysmic binaries, a magnetic 
white-dwarf mass accretes matter from its Roche-lobe filling 
main-sequence companion. Part of the energy released in the 
process appears as soft X-ray emission of the heated photo- 
sphere in the accretion spot(s) near the magnetic pole(s) of the 
white dwarf. The emitted spectrum is conventionally modeled 
as a single blackbody with a temperature ar ound 30 eV (e.g. 
iRamsav et aiTl 996: Ra msav & Cropperll2004l) . although on the- 
oretical grounds the source is not expected to be isothermal. 
Constructing a viable model requires observations that extend 
from the soft X-ray into the far-ultraviolet regime with suffi- 
cient spectral resolution and good counting statistics, prefer- 
ably from a source that exhibits a single accretion spot. A suit- 
able source is the prototype polar AM Her in its high state, 
when it is the brightest soft X-ray source in the sky. The first 
high-resolution sof t X-ray spectrum of AM Her was taken by 
Paer els et al.l (fl996) with the EUVE short-wavelength spectrom- 
eter. Because of the limited spectral coverage, however, they re- 
tained the single-blackbody assumption. Substantially improved 
high-resolution soft X-ray spectra of AM Her became avail- 
able with the Low Energy Trans mission Grating Spectrometer 
(LETGS) on board of Chandra dBurwitz et al]l2002t lTrillll2006l) 
with the High Resolution Camera for Spectroscopy (HRC-S) as 
detector. In this paper, we analyze the orbital mean LETG spec- 
trum of AM Her taken in September 2000. 

Intermittent heating and cooling determines the mean tem- 
perature distribution averaged over time and the spatial coordi- 
nates in the spot. The distribution range s from the photospheric 
temperature to several tens of eV (see iBeuermann et al.l l2008. 
for the case of AM Her). We determined the mean tempera- 
ture distribution for AM Her by unfolding the time-averaged 
Chandra LETG spectrum of AM Her. Our model employs a 
quasi-continuous distribution of blackbodies (BB) that is de- 
scribed by a single free parameter besides the common nor- 



malization. Mathematically, this model is as simple as the sin- 
gle blackbody. We fitted the LETG spectrum of AM Her by 
adding the thermal emission of the post-shock cooling flow to 
the soft X-ray comp onent. A similar model was adopted by 
dTraulsen et al.ll2010l) for the XMM Newton spectrum of the po- 
lar AITri. 



2. Observations of AM Her and data analysis 

We analyzed the Chandra observation 645 of 30 September 
2000, which was performed with the LETG and the HRC-S 
as detector and lasted 24 ks. The data cover the wavelength 
range 4-120A (0.1-3.0 keV) with a spectral resolution of 
A/t(FW HM) 0.05 A . The time-dependent spectra were ana- 
lyzed bv lBurwitz et al.1 (|2002) and only the time-averaged mean 
is considered here. The spectrum was re-binned from its original 
resolution of 0.05 A to a logarithmic wavelength scale with 100 
bins per decade, which averages over the multitude of faint emis- 
sion and absorption features but leaves the major spectral emis- 
sion lines intact (Fig.[TJ left panel, solid green curve). Absolute 
spectral fluxes were derived using the LETGS effective areas 
given in the Chandra calibration files CALDB 4.4.7 Q We re- 
analyzed, in addition the ROSAT PSPC observation rp300067 
of 12/13 April 1991 dGansickeetal.i l 19951 IRamsav et al.lll996T: 
IBeuermann et al.l 12008b . This spectrum contains 1 .3 million 
counts and has a statistical error near 0.20 keV of 0.39% per 
10-eV PSPC channel. Both data sets provide nearly full orbital 
phase coverage. On both occasions , AM Her was in its normal 
accretion mode dMazeh et al.lll986l) with visual AAVSO mag- 
nitudes at orbital maximum of V - 13.3 (Chandra) and 13.2 
(ROSAT), respectively. The PSPC spectrum was a nalyzed us- 
ing th e detector response matrix DRMPSPC_A01c dBeuermannl 
2008). 



* Based on observations with the Chandra and ROSAT satellites. 



1 The calibration files can be found at http://cxc.harvard.edu/caldb/ 
downlo ads/Release_notes/CALD B_v4.4.7,html. For a general approach 
see also Beuermann et al. (2006). 
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Energy Energy (keV) 

Fig. 1. Left: Mean Chandra LETG spectrum of September 2000 reduced to the first diffraction order (green curve) and fitted with 
the EXP+CF model (upper black solid curve, see text). Right, top: Residuals for the 1BB+CF model represented by the data/model 
flux ratio. The reduced^ 2 value refers to the energy range of 0.10-0.50 keV. Center: Same for the 2BB+CF model. The twelve 
spectral bins excluded in the fits are indicated by tick marks. Bottom: Same for the EXP+CF model displayed in the left panel. 

Table 1. Parameters of spectral fits to the mean September 2000 Chandra LETG and April 1991 ROSAT PSPC spectra of 
AM Her. The upper index '15eV refers to spectral components with blackbody temperatures kT > 15 eV. BB=blackbody model, 
CF=cooling-flow model, EXP=exponential model (Eq. 1), RL=' He-like emission lines' excluded from the fit to the Chandra LETG 
spectrum. 
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(b) ROSAT PSPC, 
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TVh is for a hydrogen ionization x H = 0.25. 2) Fit excludes the helium-like emission lines. 3) x 2 an d the number of d.o.f refer to the 0. 10-0.50 keV 
energy interval. 4) RL component included in the fit. 5) x 2 an d d.o.f. refer to 0.09-2.00 keV. 6> Average count rate-weighted blackbody temperature. 



3. Model assumptions 

Our model includes two dominant spectral components, ther- 
mal hard X-rays from the post-shock cooling flow (CF) and 
soft X-rays produced by various modes of reprocessing of 
the energy transp orted from the CF into t he atmosphere of 
the white dwarf (iKuiipers & Pringld Il982t iFrank et al.1 fl988t 
Gans icke et al.l 1 1995b iKo nig et alj l2006l) 7ln addition, a small 
component of residual emission lines (RL) exists that is not ex- 
plained by the CF model and may b e due to photoionization 
dMukai et al.ll2003t iGirish et al.l l2007). We describe the soft X- 
ray component by single and multi-temperature variants of the 
blackbody model (IBB, 2BB, 10BB, EXP), further explained in 
the next Section. 

The thermal X-ray emission of the CF is described as op- 
tically thin free-free and bound-free emission for a plasma of 
solar abundances, employing Mekal spectra. This model ne- 



glects the influence that cyclotron emission exerts on the tem- 
perature distributio n and is appropriate f or AM Her as a low- 
field polar dFischer & Beuermann! 1200 lb . Depending on the 
specific accretion rate, the shock may be free-standing above 
the atmosphere or buried in the photosphere and a variety of 
such columns may coexist in the spot at any given time. For 
matter of solar composition, the shock temperature is T s = 
32.2(M/M Q )/(tf/10 9 cm)keV, or T s = 33.5 keV for the white 
dwarf in AM Her with mass M - 0.78 M and radius R - 
7.5 x 10 8 cm (Gans icke et alf2 006). X-rays emitted from buried 
columns can escape only by passing through atmospheric ma- 
terial with column densities ranging up to Nh ^ 10 25 H- 
atomscrrr 2 (Beardmore et al. 1995; llshida etal.ll997tlChr istian 
2000). We include a partial absorber with a column density /Ve.int 
and a covering fraction / c , but disregard the reflection com- 
ponent, which becomes prominent at higher photon energies 
(Ivan Teeseling et al.ll996l) . This simplified treatment is adequate 



2 



K. Beuermann et al.: A new soft X-ray spectral model for polars 



for the limited energy range covered by the LETGS, but can not 
be expected to provide an accurate measure of the bolometric 
flux contained in the hard X-ray component (see Sect. 4.4). 

We included photoabsorption by an interstellar column den- 
sity A^h of partially ionized matter with solar abundances and 
parts of the metals condensed into dust. The ionization fractions 
of hydrogen and helium, xu and xe e , are parameters of the fit, 
with the resulting value of A^h depending on x^. Most results are 
quoted for xh=0.25 and XHe-0-40, as deriv ed for near-solar in- 
terstellar space (see Beuermann et al. I I2006L their Sect. 4.2). We 
included the NeVI absorption edg es at 78 and 85A detected in 
the EUVE spectrum of AM Her dPaerels et al.| [l996). which are 
either of interstellar or circumbinary origin. The respective opti- 
cal depths, T78 and r§5, are free parameters of the fit. 

4. Results 

The CF component reproduces the 0.6-3.0keV part of the LETG 
spectrum and accounts for the hydrogen-like emission lines of 
nitrogen, oxygen, magnesium, neon, silicon, and sulphur, but 
largely lacks the helium-like triplets and some other lines. We 
excluded, therefore, twelve of 151 spectral bins from the fit that 
contain these lines and represent the RL component (marked 
by ticks in Fig. [1] center right panel). The fit yields a covering 
fraction f c 0.70, implying that 30% of the CF component es- 
cape unattenuated and 70% pass through an internal absorber of 
Miint - 3 . 1 x 1 23 H-atoms cm" 2 . 

Interstellar absorption by a column density A^h acts on both, 
the soft and hard X-ray components. For our standard values 
xh=0.25 and XHe=0.40, the derived column densities in Table 1 
exceed the value of (3.0+ 1.5) x 10 19 H-atoms cm ' 2 derived from 
Lyq a bsorption in the FUV spectrum of AM Her (Gan sicke et al.l 
1995). The fit value of A^h depends, however, on xh and for our 
final model the Lya value is reproduced for xu = 0.68 + 0.18. 
The true value of the degree of ionization along the line of sight 
to AM Her is not well known and we can not exclude such a 
high value. We note that AM Her itself does not add much to the 
mean ionization along the path, because its Stroemgren radius 
amounts only to a few pc for a mean atomic hydrogen density 
of 0.25 hydrogen atoms cm' 2 . Our fits confirm the absorption 
edges at 78 and 85 A found by iPaerels et all (Il996l) . although 
with smaller optical depths T78 0.05 and ^ 0.11. 

Obtaining a formal x 2 °f the fit to the LETG spectrum 
requires an assessment of the errors of the binned spectral 
fluxes. This is not straightforward, because faint emission and 
absorption lines are spread ov er the entire wavelength range 
dBurwitz et al.l2002tlTriilll2006l) . Since these features are not ac- 
counted for by our model, we defined the wavelength-dependent 
'error' for the x 2 calculation by a polynomial fit to the rms devi- 
ations between data and the best model vs. wavelength, with the 
twelve spectral bins defining the RL component excluded. This 
fixes the residual xl of the best-fitting models near unity. The 
adopted error varies from 5% of the first-order flux at 0. 1 keV to 
13%at3.0keV. 

4. 1 . Modeling the soft X-ray component of AM Her 

1BB+CF model: The single-blackbody model utterly fails to 
fit the low-energy part of the LETG spectrum. Our fit yields 
an unacceptable reduced xl - 4-71 (Table 1). The residuals, 
in the form of the ratio data/model, are displayed in Fig. Q] (top 
right panel). The model involves two free parameters for the soft 
X-ray component, three parameters for the CF component, and 
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Fig. 2. Projected emitting area per eV on the white dwarf in 
AM Her and integrated blackbody flux per eV at the Earth as 
a function of the blackbody temperature kr for the 10BB+CF 
model. The solid line represents the EXP model of Eq. 1 (see 
text). 



three parameters describing the interstellar absorption (A'h, t"78 
and T85, with xu and xn e fixed). The fit parameters and the de- 
rived fluxes are summarized in Table 1 . 

2BB+CF model: Adding a second blackbody (and two free 
parameters) significantly improves the fit with best-fit tempera- 
tures of 24.7 and 41 .5 eV. The residuals of the fit are displayed in 
the center right panel of Fig. Q] The fact that the 2BB +CF model 
fits the LETG spectrum so well might suggest a dichotomy in 
the blackbody temperature of the accretion spot. On theoretical 
grounds, we consider it more plausible, however, that the result 
indicates the presence of a wider range of blackbody tempera- 
tures in the time-averaged spectrum. 

10BB+CF model: We study the temperature distribution in 
more detail, by introducing a grid of 10 fixed blackbody temper- 
atures, which covers the interval 10-60eV in 5 eV steps with 
kT n = 2.5 (2n-l)eV and n-3 to 12. The fit turns out to be stable 
if a small amount of regularization in the normalizations is in- 
troduced. We left the normalizations for nine of the ten bins free 
and prevented a loose-end excursion in the first bin by tying it to 
the trend defined by the other normalizations. Fig. [2] shows the 
best-fit projected emitting area per eV, a{T), and the correspond- 
ing unabsorbed integrated blackbody flux per eV, fbb(T). The 
latter peaks in the 15-20eV bin and substantial contributions are 
expected to exist at still lower temperatures. The observed soft 
X-ray spectral flux lacks these low-temperature contributions be- 
cause of interstellar absorption. Only a meagre 1.2% of the inte- 
grated observed flux originate from kT < 15 eV. The value of x 2 
is slightly improved for the 10BB model compared with the 2BB 
case, but xl deteriorates from 1.03 to 1.13 (seeTable 1), because 
the individual normalizations are not independent of each other. 
The errors in Fig. |2]are estimates from runs with different levels 
of regularization of the first bin. We conclude that the emitting 
area on the white dwarf decreases continuously with increasing 
temperature, although with ten free parameters the 10BB model 
is clearly overdetermined. 
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EXP+CF model: For kr > 15 eV, the projected emitting area 
a(T) of the 10BB model can be approximated by an exponential 



1000 



a{T) = a exp(-r/r ) 



(1) 



shown as the solid line in in Fi g. [2] For a distanc e of AM 
Her of 80 pc (iThorstensenl 12003; Be uermann! 120061) . we find 
a = 1.17 x 10 17 cm 2 eV~ 1 and kT = 3.46eV. The standard ver- 
sion of our EXP model utilizes the full temperature range and 
fits the soft X-ray emission of AM Her excellently with only one 
free parameter besides the normalization. Hence, the EXP model 
is mathematically as simple as the single blackbody. A fit using 
this model is displayed in Fig. [TJ (left panel and bottom right 
panel). The parameters and derived fluxes are given in Table 1 . 

The EXP model is attractive, because it provides a heuris- 
tic description of the time-averaged mean temperature distribu- 
tion in the photospheric accretion spot on the white dwarf in 
AM Her. The theory of the heating and cooling processes in the 
individual subcolu mns of the spot is not yet well developed (e.g. 
Koni g et aDl2006h and the present analysis may help to improve 
it. 

As shown in the Appendix, some key quantities of the EXP 
model can be expressed analytically. For the temperature range 
from kr = 15 eV to infinity (or ~ lOOeV) and ao and k7o as 
derived from the best fit, the projected area of the soft X-ray 
emitting spot is A±x = 4.8 x 10 15 cm 2 at of = 80pc, corresponding 
to 0.07% of the surface area of the white dwarf. For still lower 
temperatures, the spot area increases, but can no longer be deter- 
mined from the fit to the LETG spectrum alone. 

4.2. The FUV flux of AM Her 

In the EXP model, temperatures kr ^ 15 eV contribute lit- 
tle to the soft X-ray flux detected with the LETGS, but pro- 
duce a sizeable fraction of the FUV flux. In principle, the run 
of a(T) at low temperatures can be obtained from FUV obser- 
vations, subject, however, to substantial observational and the- 
oretical uncertainties. Observationally, the FUV flux is domi- 
nated by emission from the accretion stream, with the Balmer 
(and presumably the Lyman) jump strongly in emission. The 
spot contribution, on the other hand, is probably character- 
ized by rapid variability and the lack of strong line emission. 
Whether it shows t he Lyman jump in emission or absorption 
remains uncertain (Gansicke et alj Il995| IT998I; iGreelev et al.l 
1999; iHutchings et alJl2002t iKonig et aljhoOrjh . For simplicity, 
we retain the multi-temperature blackbody model. 

The mean orbital spectral flux of AM Her at 930A was 
measured with the shuttle- based Hopkin s Ultra violet Telescope 
(HUT) in March 1995 (IGreelev et al.l [l999b and with the 
Far Ultraviolet Spectroscopic Explorer (FU SE) in June 2000 
dHutchings et all 120021: iGansicke et al.l 120061) as / 930 - 3.5 x 
10~ 13 and 2.6 x 10" 13 ergcm-V'A" 1 , when the mean AAVSO 
visual magnitudes of AM H er were V - 13.2 and 13.3, respec- 
tively. IGreelev et al.l (1 19991) reported the presence of a rapidly 
fluctuating 'flare minus non-flare' component (FNF) with a 
mean-orbital 930 A spectral flux of 1.3 x 10~ 13 ergcm _2 s _1 A _1 . 
Corrected for reddening, the quoted fluxes become /930 - 
4.0 x lO^ergcnrV'A- 1 (HUT), 3.0 x lO" 13 ergcnrVA- 1 
(FUSE), and 1.5 x 10~ 13 ergcnrVA- 1 (HUT, FNF). We com- 
pare these observed fluxes with the predictions of the spot mod- 
els presented in the last Section. We do not correct for the or- 
bital modulations of the individual components, because the ob- 
servations are not simultaneous, the individual high states have 



> 



C 
3 
O 

U 



100 



1 1 I 



~i — 1 — 1 — r~| 

AM Her 

PSPC Apr 1991 
EXP+CF+RL 




0.2 0.5 
Energy (keV) 

Fig. 3. Top: 2BB, LIN, and CF components of the September 
2000 Chandra LETG spectrum of AM Her fitted to the mean 
April 1991 PSPC spectrum. Bottom: Residuals of the fit based 
on the purely statistical errors of the PSPC spectrum. 

different flux levels, and the blackbody assumption adds to the 
uncertainties. 

The EXP model yields a 930 A spectral flux / 930 = 8. Ox 10" 13 
erg cm _2 s _1 A -1 if integrated over all temperatures and f^ v - 
1 .35 x 10~ 13 erg cm~ 2 s -1 A -1 for kr > 15 eV. The former exceeds 
the reddening-corrected total HUT and FUSE 930 A fluxes by 
a moderate factor of two, while the latter agrees approximately 
with the lower HUT FNF flux. Taken at face value, these num- 
bers suggest that the EXP model overestimates the emitting 
area at kr < 15 eV. Reducing a(T) by a factor of four at 
5 eV (dashed line in Fig. |2]i yields an acceptable /930 = 4.06 x 
10~ 13 erg crn-VA" 1 , while f^ v = 1.33x10-° erg crrrV 1 A" 1 
stays practically unchanged. These are the numbers included in 
Table 1 . Since the inherent uncertainties in this comparison may 
easily amount to a factor of two, we conclude, that the EXP 
model or its variant with a reduced emitting area for kr < 15 eV 
provide an adequate description of the FUV flux expected from 
the entire hot polar cap of AM Her. The fraction of the flux that 
arises from kr > 15 eV, agrees with Greeley's FNF component, 
suggesting that it represents the Rayleigh- Jeans part of the spec- 
tral flux emitted by the proper accretion spot. For comparison, 
the IBB model fails to predict the FUV flux by more than an 
order of magnitude (Table 1). 

4.3. Fitting the ROSAT PSPC spectrum 

Finally, we show that the EXP+CF+RL model fits the April 
1991 ROSAT PSPC spectrum, too. We allowed for small cor- 
rec tions to the nom inal calib ration of the PSP C as prescribed 
by iPrieto et aD (1 19961) and iBeuermannl (120081) . There are four 
free parameters of the fit, two for the EXP component and two 
for the normalizations of the numerically provided RL and to- 
tal CF components taken from the fit to the Chandra LETG 
spectrum. The excellent fit is shown in Fig. [3] and the derived 
parameters are listed in Table 1 . For the PSPC fit, the normal- 
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Table 2. Bright-phase high-state bolometric energy fluxes of 
AM Her 



Component 



(erg cm 2 s ') 



r\n L Ref. 

(ergs" 1 ) 



Soft X, BB kr> 15 eV 


5.0 


X 


10" 


9 


In 


1.9 


X 


10" 


1 


Hard thermal X, CF 


3.1 


X 


io- 


10 


3.37T 


2.0 


X 


10 32 


2 




6.4 


X 


io- 


lo 


3.3;r 


4.0 


X 


10 32 


3 


Cyclotron 


1.1 


X 


io- 


lo 


In 


4.2 


X 


10 31 


4 


Accretion stream 


1.7 


X 


io- 


id 


4n 


1.3 


X 


10 32 


4 



References: (1) LETG EXP fit, this work; (2) LETG CF fit, this work, 
AWmai = 3 x 10 22 H-atoms cm- 2 ; (3) RXTE, Christian (2000), Mtintemai 
up to 10 25 H-atoms cm -2 ; (4) Optical/IR, Bailey et al. (1988), Gansicke 
et al. (1995). 

ization of the EXP component is almost identical to that of the 
LETG, while the CF-component is 24% weaker and the RL com- 
ponent 27% stronger, differences that are easily accounted fo r 
by time variability of the source (IStella et al.ll 19861; iTrilll 120061) . 
An acceptable fit to the PSPC spectrum is also obtained for the 
1BB+CF+RL model (Table 1), demonstrating that the higher re- 
solved LETG spectrum is needed to define a superior model for 
the soft X-ray emission of AM Her. 



4.4. Accretion luminosity and accretion rate 

Soft and hard X-ray emission are the dominant accretion- 
induced radiation components in AM Her. Table 2 summarizes 
their integrated fluxes, which refer to the bright phase (magnetic 
phases mag n = 0.2 - 1.0) and exceed the mean orbital fluxes in 
Table 1 by about 22%. The fluxes F are converted to luminosi- 
ties L = rjnd 2 F, assuming the geometry factors 77 given in the 
table and a distance d — 80 pc (Thorstensen 2003t lBeuermannl 
2006). The ultraviolet flux from the heated polar cap of the white 
dwarf is not included, because it is created by the downward 
hard X-ray and cyclotron fluxes dGansicke et al.l 1 19951 Il998t 
iKonig etalJ l2006) and is already accounted for by the choice 
of the geometry factors. The fact or 33n for the hard X-ray s ac- 
counts for the reflection albedo (Ivan Teeseling et aTI [l996). We 
adopt here Christian's (2000) integrated hard X-ray flux. The to- 
tal accretion-induced luminosity then is L acc = 2.5 x 10 33 erg s _1 
and the implied accretion rate for a white dwarf with a mass 
of 0.78 M and a radius of 7.5 x 10 8 cm dGansicke et al.i r2006) 
is M = 2.9 x lQ-^Mnyr- 1 . Thi s result is similar to the PSPC 
result of lBeuer mann et alj (120081) . but is placed here on more se- 
cure grounds, being based on the much better resolved Chandra 
LETG spectrum. The luminosity emitted as soft X-rays exceeds 
the sum of the hard X-ray and the cyclotron luminosities by a 
factor L S oftx/(ihardx + ^cyc) - 4.3 + 2.0. The conservatively esti- 
mated error accounts for a ~ 30% uncertainty in the hard X-ray 
flux, the remaining error in the bolometric soft X-ray flux, and 
variations between the individual high states in which the mea- 
surements were taken. A moderate dominance of soft X-rays is 
clearly p resent and is easily explained by the concep t of 'buried 
shocks' (Kuijp ers & Pringldl 19821; [Frank et al.lfl988b . 



5. Conclusions 

We have presented a new model for the soft X-ray emission of 
polars, which is superior and physically more plausible than the 
frequently employed single-blackbody model. Our model pro- 
vides an excellent fit to the high-resolution Chandra LETG spec- 
trum of AM Her and fits the ROSAT PSPC spectrum equally 



well. Without further assumptions, the model provides also an 
explanation of the ra pidly varying FUV spectral flux of AM Her 
(iGreelev et al.ll 19991) . which may represent the Rayleigh Jeans 
tail of the soft X-ray component. On theoretical grounds, the 
intermittent heating and subsequent cooling of individual seg- 
ments of the accretion spot must lead to a spatially and tem- 
porally averaged temperature distribution in the spot, which ex- 
tends from several tems of eV down to the photospheric temper- 
ature. Our model provides a first approximation to this temper- 
ature distribution, using a single-parameter fit to the soft X-ray 
spectral energy distribution. The new model is mathematically as 
simple as the isothermal single-blackbody model and physically 
more plausible. 
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Appendix A: Analytical expressions 

Let the emitting area be a(T) = a^exp(-T/To) and the emitted 
spectrum at temperature T be a blackbody. For a temperature 
interval of T\ to infinity and x = T/T^, the projected area of the 
spot is 



I 



A j. = Q.d - flo To I e v dx = an ?o e 



(A.l) 



with O the solid angle subtended by the spot and d the distance. 
With C — oo crT^/ind 2 ), where <x is the Stefan-Boltzmann con- 
stant, the integrated unabsorbed blackbody flux at the Earth is 

o 

F bh = C I *V* dx = C(x\ +5 x] + 12 xf +24 x x +24) t~ x \ (A.2) 



The mean spot temperature weighted by the unabsorbed black- 
body flux is 



(T hb ) = T 



5 + 



x\ + 5 x 3 + 12 x\ + 24 xi + 24 ) 



(A3) 



which assumes its minimal value (7bb) = 5 To if me exponential 
extends to zero temperature. 
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